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ABSTRACT 

Recent observations show a strong correlation between stellar mass and accretion rate in young 
stellar and sub-stellar objects, with the scaling Mace oc holding over more than four orders of 
magnitude in accretion rate. We explore the consequences of this correlation in the context of disk 
evolution models. We note that such a correlation is not expected to arise from variations in disk 
angular momentum transport efficiency with stellar mass, and suggest that it may reflect a systematic 
trend in disk initial conditions. In this case we find that brown dwarf disks initially have rather larger 
radii than those around more massive objects. By considering disk evolution, and invoking a simple 
parametrization for a shut-off in accretion at the end of the disk lifetime, we show that such models 
predict that the scatter in the stellar mass-accretion rate relationship should increase with increasing 
stellar mass, in rough agreement with current observations. 

Subject headings: Accretion, accretion disks — stars: pre-main sequence — stars: low-mass, brown 
dwarfs — planetary systems: protoplanetary disks 
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1. INTRODUCTION 

It is well established that young stars of around 
solar mass, such as T Tauri stars (henceforth TTs), 
are surrounded by disks. Reprocessing of stellar radi- 
ation by these disks causes excess emission above the 
stellar photosphere in the infrared, and the accretion 
of material on to the stellar surface produces both 
line emission and excess blue continuum emission. 
This permits an observational measurement of the 
instantaneous ste llar accretion ra te for large samples 
of TTs (e.s iCuUbring et aLllTOQl . More recently, it 
has become clear that young brown dwarfs (henceforth 
BDs) also posses circumstellar dis ks, with ob servations 
detecting infrared excesses (e.g. INatta fc Te sti 2001; 
Javawardhana et al. '2003b ^, acc re tion si g natures (e.g. 
Javawardhana. Mo hantv & Basril l2003al: INatta et al.l 
20041 iMohantv. J avaward hana fc Basrf I2005 D and 
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millimeter continuum emission (e.g. iKlein 

If BD formation is merely a scaled-down version of 
star formation, one might guess that the disk mass 
and stellar accretion rate Mace should scale roughly 
proportional to the (sub-)stellar mass M*. In fact, 
the accretion rate on to the "stellar" surface shows a 
striking correlation wi t h the square of th e stellar mass 
jMuzerolle et al. 2003; Natta et al. 2004; 'Calve t et all 
12004: Mu zcroUc ct al. 2005; Mo hantv et al.. .20051 The 
correlation shows a large scatter, but holds over more 
than 2 orders of magnitude in mass and 4 orders of 
magnitude in accretion rate. We expect a large scat- 
ter, as disk accretion r ates decrease as disks evolve (e.g. 
IHartmann et al.l I1998|) , but the physical origin of the 
Mace oc M^ relationship is not clear. Here we suggest 
that the most conservative interpretation is that the cor- 
relation reflects the initial conditions established when 
the disk formed, followed by subsequent viscous evolution 
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of the disk. We explore the observational consequences 
of this model. 

2. MODEL 

The equation for the surface density evolution, S(i?, t), 
of a viscous accretion disk is 



as 

dt 



3__d_ 

RdR 



dR 



(1) 



where R is cylindrical radius, t is time, and v{R,t) is 
the kinematic viscosity of the disk |Pringlc 1981). If the 
viscosity is independent of time and can be expressed 
as a power-law function of radius v oc , then Equa- 
tion^permits exact analy tic solutions. Here we co nsider 
the similarity solution of'Hartmann et all l)1998fl . after 
^fiidcn-Bcll & Pringlc (1974). In this case the solution 
for the surface density of the evolving disk is 



Md(0)(2-7) -(5/2-T) 

I](i?, t) = ^_ n2„., T exp 



27ri?2r7 



(2) 



where Md{0) is the initial disk mass. The dimensionless 
radius r = R/Rq, where Rq is a scale radius which sets 
the initial disk size. The dimensionless time t = t/t,y + 1, 
where the viscous scaling time is given by 



3(2-7)2^^0 
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Here Vq is the value of the viscosity at radius Rq. The 
accretion rate on to the star is therefore 



Mace = 



MdjO) 
2(2 -7)t,' 



(4) 



We suggest that the time evolution implied by this equa- 
tion is responsible for some of the observed scatter in the 
mass-accretion rate distribution, and assume that the ob- 
served Mace OC Ml relationship holds for the initial ac- 
cretion rates. Observations of BD disks have shown that 
the disk- to-star mass ratio is comparable to that found 
for TTs ijKlein et al.ll2003l) . so we assume that the initial 
disk mass scales linearly with the stellar mass. In order 
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to reproduce the observed correlation, that Ma 
we therefore require that the 



1 

Ml 



(5) 



For TTs, t he viscous time-scale ty is typically ~ 10''- 
lO^yr fe.g. IHartmann et alJUQQ^ . and therefore we ex- 
pect that the viscous scaling time for BDs be of order 
lO^yr. 

However we no te that while current observations (e.g. 
iKlein et 3111200.111 suggest a linear relationship between 
disk and stellar mass the data do not yet rule out a 
steeper relationship. A much steeper scaling of cx 

could also explain the observed Mace oc M^ relationship 
without requiring the change in with stellar mass seen 
in Equation 13 However in this case BD disks are less 
massive (relative to the central object) then their TT 
counterparts, with few BDs having disks more massive 
than ~ IMj. This seems unlikely, but current observa- 
tions do not yet rule it out. 

We note that it is also possible to satisfy the observed 
relationship if t^, increases with increasing stellar mass, 
so that Mace cx M^ only at late times (i.e. t ty). We 
reject this scenario for three reasons. Firstly, it requires 
extremely short viscous time-scales for BD disks {ty < 
lO^yr), which seem unlikely. Also, this scenario requires 
that the initial accretion rate be higher for BDs than 
TTs, which again seems unlikely. Lastly, in this case 
most of the initial disc mass is accreted very rapidly. By 
t IMyr BD disk masses are only a few percent of a 
Jupiter mass, which is rather less massive than observed 
l|klein et al...2003,l . 

2.1. Viscosity 

We now consider the consequences of this relationship 
in terms of a simple alpha -model for the disk viscosity 
ijShakura fc Sunvaevlll 97,111 . Theoretical studies of angu- 
lar momentum transport in disks suggest that magne- 
tohydrodynamic (MHD) turbulence can produce values 
of a of order 0.01 fe randcnburg et al. 1995; Stone et a^ 
119961 iBalbus fc HawlevI 119981) . ' In a Keplerian disk we 
can re-express this form for the viscosity in terms of the 
disk midplane temperature T and the angular frequency 

T TR^f^ 

V oc OL— cx OL 777^ , (Oj 
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If we substitute the viscosity at the scale radius i?o into 
Equation 13 we find that the viscous scaling time of the 
disk satisfies 



ty oc 



rTmV^ 

aT{Ro) 



(7) 



MHD simulations suggest that a does not vary signif- 
icantly with M». Also, constant aspect ratio H/R de- 
mands that the disk midplane temperature scales as 
T oc M^/R. Therefore in the case of constant aspect 
ratio H/R we can solve for the dependence of Rq with 

M, exactly to find i?o oc MZ^'^ . In real TT and BD 
disks, however, reprocessing of stellar radiation is the 
dominant mode of disk heati ng and typically results in a 
"flaring" disk structure (e .g. IKenvon fc Hartmannlll98'li 
IChiang fc Goldreichlll997fl . In this case the temperature 



scales as T oc i?^-'/^, which in turn gives Rq oc A/* ^ . 
The exact form of the disk midplane temperature is 
therefore somewhat model-dependent, but it is clear that 
in order to satisfy Equation|31we require that Rq decrease 
with increasing stellar mass. This implies that BD disks 
are initially larger than those around TTs, with typical 
scale radii Rq for BD disks of around 50-lOOAU. This 
is not altogether unexpected, as if the material in a col- 
lapsing cloud has a roughly constant specific angular mo- 
mentum we expect smaller disks around more massive 
objects. We note that the disks expand in size on their 
viscous time-scale, and therefore TTs disks will rapidly 
expand to be comparable in size to BD disks. Addition- 
ally, externa l factors, such as tidal truncation by a binary 
companion l)Paczvhskilll977l) . can strongly influence the 
size of an evolving disk. It is therefore unlikely that this 
predicted scaling of the initial disk size with stellar mass 
will have a significant observational signature. 

2.2. Cessation of accretion 

Observations of TTs show that the the cessation of 
accretion occurs concurrently with (inner) disk clear- 
ing (e.g. iHartigan et alJ 1799(11 . and the data are con- 
sistent with individual objects undergoing a rapid 
"shut-off" in acc retion once the accretion rate fa lls to 
some low level () Armitage. Clarke fc Pallal 12003(1 . It 
has been suggested that photoevaporation by the cen- 
tral star is responsible for this shut-off in accretion 
(fClarke, Gendrin fc Sotomavor 2001), with accretion ter- 
minating once the disk accretion rate falls to a level com- 
parable to the expected photoevaporative disk wind rate. 
The integrated mass-loss rate from a photoevaporative 
disk wind is given by 



1.6xlO-i"Moyr-i (^^^ 
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(8) 

where 4> is the rate at which the star emits ioniz- 
ing photons (jllo llcnbach et al. 1994; Font et al. 2004). 
lAlexander. Clarke fc Pringld lj2005ll show that chromo- 
spheric activity in TTs can provide an ionizing flux 
of order lO^'^photons s~^ or higher. The photospheric 
emission from both TTs and BDs is powered by en- 
ergy released from gravitational collapse, but little is 
known about the magnetic activity of BDs. However 
recent observations have found that the ratio of the X- 
ray:bolometric lu minosity is approxima tely the same for 
BDs as for TTs l|Preibisch et al.ll2005|l . so we make the 
crude assumption that 4> scales linearly with mass. Thus 
the predicted photoevaporative wind rate for a O.OIM© 
BD is of order IQ-i^Mgyr-i. 

2.3. Model parameters 

We have constructed disk evolution models in which 
the initial disk mass is given by Mdifi) = O.OIM*, and 
the viscous scaling time of the disk is given by 



t. 



5 X lO^yr 



IM^ 



(9) 



We note, however, that some scatter in both Md{Q) and 
ty is likely. In our fiducial model we adopt a viscosity in- 
dex 7 = 1, and we have evaluated models across the mass 
range = 0.01-IOMq. 7 = 1 is consistent with a T cx 
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Fig. 1. — Evolution of accretion rate in viscous disk models 
with varying stellar mass. The models have initial disk masses of 
O.OlAf* and viscosity index 7 = 1, and the viscous scaling time 
for the fiducial IMq model is 5 X lO^yr. The heavy solid line 
shows how the accretion rate in the fiducial model is shut-off by 
a photoevaporative wind. The heavy dashed line shows how the 
evolution of the accretion rate changes if a different viscosity index 
(7 = 3/2) is adopted (with the same initial accretion rate). 

R^^l'^ flared disc model of t he type discussed above (e.; 
IKenvon fc HartmannI IToSTI IChiang fc GoldreichI [1997 

but 7 is not especially well-constrained by current ob - 
servations (see e.g. discussion in Hartma nn et a 
FigQ] shows how the accretion rate in the models evolves 
with time for different stellar masses and viscosity laws: 
we see that the accretion history of these models is not 
strongly dependent on the choice of 7 over the timescales 
in which we are interested. Also shown in FigQ is the 
evolution of a model which incorporates disk photoevapo- 
ration (using the models of Alexander. Clarke & Pringle 
I2OOI . In this case the accretion rate is unchanged over 
most of the evolution, but accretion is rapidly shut-off at 
the end of the disk lifetime. 

We have evaluated the above model set but assumed 
that accretion is shut-off once the accretion rate falls 
below some critical value. We define this value as in 
Equation 121 and assume that the ionizing flux scales as 
$ = lO^'^s"^ (M^/IMq). FigElshows how accretion rate 
varies as a function of stellar mass in these models, and 
also how the accretion rates change as the disks evolve. 
For comparison we include observed data points from a 
number of recent studies. Our simple model reproduces 
both the slope and scatter in the observed distribution 
well, from the BD regime up to masses of several Mq . 

The assumptions in our model lead to some notable ob- 
servable consequences. The longer viscous scaling time 
for disks around lower mass objects means that the evo- 
lution of these disks is slower (see Fig^. However the 
assumed scaling of $ with stellar mass means the rate at 
which accretion is shut-off decreases only linearly with 
stellar mass, and so BD disks evolve for fewer viscous 
times than TT disks before accretion ceases (see FigEJ. 
As a result we expect a smaller spread in observed ac- 
cretion rates as stellar mass decreases, with BD accre- 
tion rates showing less scatter than those for TTs. This 
is in tentative agreement with the data shown in Fig[21 
where the detected accretion rates show a (2-cr) scatter 
of il.Odex in the mass range O.O2-O.2M0, but a larger 
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Fig. 2. — Accretion rate plotted as a function of stellar mass. 
Data points are taken from[Calvet e~a l. (2004, filled circles with 
error bars) , liN atta et aJJ 1200^ cross es" for detections, arrows for 
upper limits), Muzcro Ue et al.l <200l filled squares for detections, 
filled triangles for upper limits), and IMohantv et al.l 11200.1 open 
circles for detections, open triangles for upper limits). The dia- 
monds show the compilation of previous data given in Table 4 of 
Mohanty et al. (2005). The solid lines are isochrones for our disk 
model, plotted (from top to bottom) at t = 0, 0.1, 0.5, 1.0, 2.0 and 
S.OMyr. The shaded area indicates the range of predicted accretion 
rates when accretion is shut-off by photoevaporation at late times. 

scatter of ±1.4dex in the range O.4-4.OM0. The model 
also predicts that few, if any, objects should be detected 
with accretion rates significantly below 10^ ^^Mq. How- 
ever we note that this value is near to the sensitivity 
limit of current observations, as seen by the locus of the 
many non-detections (upper limits) in Fig[21 In this sce- 
nario we expect these objects to have little or no ongo- 
ing accretion, essentially appearing as "weak-lined" BDs. 
Lastly, in our model both TT and BD disks have typical 
lifetimes of order a few Myr, with no significant differ- 
ence in disk lifetimes across the mass spectrum. Cur- 
rent observations tentatively agree with this prediction 
( Mohanty et al . 2005.) , but more data are needed in order 
to make a more detailed comparison. 

3. DISCUSSION AND SUMMARY 
In the preceding Section we have outlined the observa- 
tional consequences if the M^-Macc relation for T Tauri 
stars and brown dwarfs is set by disk initial conditions 
followed by viscous evolution. Alternatively, the ob- 
served relation might reflect systematic variations in the 
structure or feeding of disks. One possibility is that some 
aspect of disk angular momentum transport introduces a 
bottleneck into the accretion flow, restricting the stellar 
accretion rate to the observed Mace oc Ml independent 
of anything happening at larger radii. [Mohanty et a^ 
I^D5) suggest that this could be caused by suppres- 
sion of the magnetorotational instability (MRI) due to 
a reduced ionization fraction in the disks around cooler, 
lower-mass stars, a scenario is similar to the layered disk 
model introduced bv ICammic (1996). Except for very 
close to the star, accretion occurs via an externally ion- 
ized surface layer of column S^, at a rate Mace oc a^S^. 
In the simplest versions of such models the dominant ion- 
ization source for the surface layer is cosmic rays, with 
the consequence that both ~ 10^ gcm~^, and the re- 
sulting stellar Mace, are constants independent of stellar 
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mass. A variant of such a model might yield the observed 
M^-Macc relation, but it does not seem promising. 

A more radical possibility is that the disk acts only 
as a short-term buffer between the star and gas con- 
tinually captured from the surrounding envir onment via 
Bondi-Hoyle accretion ijBondi fc HovlellT94^ . This im- 
mediately predicts the observed Mace oc scaling 
ijPadoan et al.l I2'00,'tD . For such a model to work, the 
captured gas must circularize at a small enough radius 
that the viscous timescale is short compared to the ~Myr 
disk lifetime, and we would expect the accretion rate 
to vary with the cir cumstellar environment. However 
iMohantv et al.l l)2005j) claim that many of the objects for 
which this correlation is seen to hold show similar ac- 
cretion rates despite a broad dispersion in their environ- 
ments. Also, for Bondi-Hoyle accretion in a turbulent or 
clumpy medium, we expect the scatter in accretion rate 
to increase with decreasing stellar mass, as the smaller 
accretion radius of lower mass stars samples a smaller 
volume in which random fluctuations in density would 
be greater. 

In this Letter, we have examined the consequences that 
follow if the observed Mace oc M^ relation for brown 
dwarfs and low mass stars reflects a combination of vary- 
ing disk initial conditions with stellar mass and viscous 
disk evolution. A relation of this form requires that the 
"initial" disk size increase with decreasing stellar mass. 
A larger initial disk size for brown dwarfs is not unreason- 
able if the specific angular momentum of the last mate- 
rial to be accreted is independent of stellar mass, though 
the specific form of the required scaling does not follow 
from any such simple consideration. Consequently we 
predict that disks around isolated BDs should be rather 



large, typically 50-lOOAU in radius. Future observations 
of disk sizes using millimeter interferometry should prove 
invaluable in testing this prediction. 

If, as we advocate, the observed Mace scaling is 

the result of disk initial conditions, we have shown that 
both the scatter in the disk accretion rate for the popu- 
lation as a whole, and the change in mean accretion rate 
with stellar age, ought to be smaller for brown dwarfs 
than for more massive pre-main-sequence stars (such as 
TTs). This provides an observational test for distinguish- 
ing between this model and alternative s in which the 
relation arises from either disk physics l lMohant v et alJ 
12005(1 or Bondi-Hoyle accretion fPad oan et al.ll2005j) . 

Lastly, following on from Section [2 we note that cur- 
rent data do not yet rule out a steeper M^,-Md rela- 
tionship than the linear scaling we have assumed. A 
much steeper scaling of Md oc M^ could also explain the 
observed accretion rate-stellar mass relationship with- 
out requiring the (spatially) large BD disks predicted by 
our model. However in this case BD disks should have 
rather low masses, with few BD disks more massive than 
^ IMj. Thus in order to for viscous disk accretion to re- 
produce the observed Mace oc relationship BD disks 
must either be spatially large, as we suggest, or have low 
masses. Future observations of the sizes and masses of 
BD disks will be invaluable in resolving this degeneracy. 
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